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1  Statement  of  the  problem  studied 


This  proposal  concerns  modeling,  analysis,  and  computational  simulations  of 
some  non-linear  interface  problems  and  applications  that,  we  believe,  are  of 
interests  to  DOD/ARO  missions  especially  in  the  fields  of  automatic  target 
recognition  (ATR),  computational  geometry,  numerical  methods  for  stochas¬ 
tic  differential  equations,  computational  fluid  mechanics,  non-linear  dynam¬ 
ics,  imaging  processing,  and  inverse  scattering. 

Most  of  our  previous  work  has  been  focused  on  linear  interface  problems 
with  discontinuities.  We  have  been  focusing  on  non-linear  interface  problems 
and  applications  that  are  of  interests  of  DOD/ARO  during  the  course  of  the 
project. 


2  Summary  of  the  most  important  results 

We  have  successfully  developed  the  immersed  interface  method  for  several 
non-linear  interface  problems  including  Poisson-Boltzmann  equations  on  ir¬ 
regular  domains,  some  variational  problems,  augmented  approach  for  inter¬ 
face  problems  and  problems  on  irregular  domains. 

2.1  Immersed  Interface  Method  for  a  non-linear  mod¬ 
eling  of  MR-fluid 

The  MR-fluid  model  is 


V  •  (/3(x,  |V«|)V«)  =  /.  (1) 

The  /3(x,  JVu.J)  represents  a  nonlinear  constitutive  law  and  depends  on  the 
local  material  property  and  composition.  For  example,  in  the  application 
of  magneto-rhological  (MR)  fluid,  say  there  are  iron  particles  within  a  fluid, 
using  a  nonlinear  saturation  model,  u  is  the  magnetic  potential  and  the 
empirical  constitutive  law  (Frohlich-Kennelly  relation)  of  the  magnetic  per- 
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meability  of  the  inclusion  is  given  by  (in  non-dimensionized  form) 


/3(x,  H ) 


< 


1  + 


(ji  -  1)  Bsat 
(n  -  1)  H  + 


if  x  is  within  the  iron, 
if  x  is  within  fluid, 


(2) 


where  Bsat  is  a  maximum  attainable  magnetization  of  the  inclusion,  j3+  is  a 
constant,  and  jU  is  the  permeability  of  iron. 

Across  the  interface,  we  have  the  natural  jump  conditions  if  there  is  no 
source  or  sink  present,  that  is 


[u]  =  0,  \/3un\  =  0.  (3) 

where  un  =  =  Vit  •  n  is  the  normal  derivative  of  the  solution  it,  n  is  the 

unit  normal  derivative. 

For  different  problems,  the  non-linear  constitutive  law  has  different  forms. 
It  is  impractical  to  discuss  them  all  in  a  single  paper.  Therefore,  in  this  paper, 
we  will  explain  our  method  particularly  for  the  problem  (l)-(3). 

Finite  element  methods  based  on  body- fitted  grid  and  the  boundary  in¬ 
tegral  integral  method  has  been  developed  for  the  linear  model.  Because  the 
coefficient  of  (1)  depends  on  the  solution,  it  is  difficult  to  use  a  boundary 
integral  or  integral  equation  method  even  for  the  linearized  problem. 

We  have  developed  a  finite  difference  method  to  solve  the  non-linear  prob¬ 
lem  based  on  a  modified  substitution  method  in  [5].  The  linearized  (l)-(3) 
problem  is  solved  using  the  black  solver  of  the  maximum  principle  preserving 
scheme.  All  the  discretizations  are  done  with  second  order  method.  Global 
second  order  accuracy  in  the  maximum  norm  is  obtained  for  a  problem  of 
which  we  know  the  exact  solution  although  it  is  almost  impossible  to  prove 
it  theoretically. 


2.2  Theoretical  and  numerical  analysis  of  the  weighted 
minimal  surface  problem 

The  weighted  minimal  surface  problem  in  piecewise  smooth  media.  The  so¬ 
lution  to  the  weighted  minimal  surface  problem  is  continuous  but  the  deriva- 
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tives  have  a  jump  across  the  interface  where  the  medium  property  is  dis¬ 
continuous.  The  jump  condition  of  the  derivatives  has  been  derived  gen¬ 
eralized  the  Snell’s  law  in  geometric  optics  to  weighted  minimal  surfaces  of 
co-dimension  one  in  any  number  of  dimensions.  We  have  developed  a  numeri¬ 
cal  method  based  on  the  gradient  flow  and  the  maximum  principal  preserving 
immersed  interface  method  to  solve  this  nonlinear  elliptic  problem  with  jump 
conditions  in  [12]. 

2.3  A  finite  difference  method  and  analysis  for  2D  non¬ 
linear  Poisson-Boltzmann  equations 

A  fast  finite  difference  method  based  on  the  monotone  iterative  method  and 
the  fast  Poisson  solver  for  a  2D  nonlinear  Poisson-Boltzmann  equation  is 
proposed  and  analyzed.  Each  iteration  of  the  monotone  method  involves  the 
solution  of  a  linear  equation  in  an  exterior  domain  with  an  arbitrary  interior 
boundary.  A  fast  immersed  interface  method  for  Helmholtz  equations  on 
exterior  irregular  domains  is  used  to  solve  the  linear  equation.  The  monotone 
iterative  method  leads  to  a  sequence  which  converges  monotonically  from 
either  above  or  below  to  a  unique  solution  of  the  problem.  This  monotone 
convergence  guarantees  the  existence  and  uniqueness  of  a  solution  as  well  as 
the  convergence  of  the  finite  difference  solution  to  the  continuous  solution. 
We  also  compare  our  results  with  available  data  in  the  literature  to  validate 
the  numerical  method.  Our  method  is  efficient  in  terms  of  accuracy,  speed, 
and  flexibility  in  dealing  with  the  geometry  of  the  domain. 

2.4  Augmented  methods  for  solving  PDEs  with  dis¬ 
continuities  and  problems  on  irregular  domains 

There  are  at  least  two  reasons  to  use  augmented  strategies  [7].  The  first 
one  is  to  get  faster  algorithms,  particularly,  to  take  advantages  of  existing 
fast  solvers.  The  second  reason  is  that,  for  some  interface  problems,  an 
augmented  approach  may  be  the  only  way  to  derive  an  accurate  algorithm. 
Using  an  augmented  approach,  one  or  several  quantities  of  co-dimension  one 
are  introduced.  The  GMRES  iterative  method  is  often  used  to  solve  the 
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augmented  variable (s)  that  are  only  defined  along  the  interface  or  the  irreg¬ 
ular  boundary.  Our  recent  work  supported  by  the  ARO  grant  including  the 
augmented  method  for  pressure  boundary  condition  of  the  Stokes  flow  with 
a  non-slip  boundary  condition  [10],  with  discontinuous  viscosity  [9,  15],  the 
bi-harmonic  equations  [1]. 

2.5  Other  contributions 

Please  see  the  list  of  publications  of  and  other  reported  information  in  the 
on-line  form  of  the  final  report. 

Our  research  results  are  summarized  in  the  recent  published  book  [8]  in 
the  SIAM  Frontier  Book  Series.  We  have  developed  high  order  immersed 
interface  method  (IIM)  in  [4],  The  immersed  interface  method  has  been 
successfully  applied  to  variety  linear  and  non-linear  problems.  It  has  been 
applied  to  a  elasticity  systems  in  [14,  17,  3];  Stokes  and  Navier  Stokes  equa¬ 
tions  with  interfaces  and  singular  sources  in  [13,  16].  Other  theoretical  and 
numerical  analysis  can  be  found  in  [11,  6,  2,  10]. 

In  a  summary,  we  are  able  to  accomplish  almost  all  the  projective  in  our 
original  proposal  and  more.  Furthermore,  the  ARO  grant  really  helped  us  to 
get  our  research  projects  going  and  move  to  new  research  areas.  We  really 
appreciate  ARO  for  the  support. 
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